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Abstract 

We argue that a primordial black hole is a natural and unique candidate for all 
dark matter. We show that, in a smooth- hybrid new double inflation model, a right 
amount of the primordial black holes, with a sharply-defined mass, can be produced 
at the end of the smooth-hybrid regime, through preheating. We first consider 
masses < 10~ 7 Mq which are allowed by all the previous constraints. We next 
discuss much heavier mass 1O 5 M0 hinted at by entropy, and galactic size evolution, 
arguments. Effects on the running of the scalar spectral index are computed. 



1 Introduction 



The presence of dark matter (DM) has been firmly established by a host of observations, 
and its abundance was measured by the WMAP with an unprecedented precision: [TJ 

tt DM h 2 = 0.1131 ±0.0034. (1) 

However it is not known yet what DM is made of, and the question remains a big mystery 
in modern cosmology as well as particle physics. 

It is often claimed that there is no DM candidate in the framework of the standard 
model (SM), assuming that DM is made of elusive particles which have evaded all con- 
ventional DM searches. The DM particle must be electrically neutral, long-lived, and 
cold, but no such particle exists in the SM. Thus we need to postulate a theory beyond 
SM and introduce a new degree of freedom, which is usually made stable by imposing an 
additional discrete symmetry. The introduction of such a discrete symmetry may be mo- 
tivated by other phenomenological reason. For instance, in the supersymmetric standard 
model (SSM), it is customary to introduce an R-parity in order to forbid dangerous op- 
erators which would give rise to too fast proton decay. Once the R-parity is imposed, the 
lightest supersymmetric particle (LSP) becomes stable, and therefore a DM candidate. 
On the other hand, there is an argument that the R-parity violation may be a common 
phenomenon in the string landscape [2]. If so, the dangerous operators must be absent 
due to some other reason(s) and the lifetime of LSP in the SSM may be too short to 
account for the DM. 

If the DM is made of a weakly interacting massive particle (WIMP), we may be able to 
observe collider, direct and indirect DM signatures; the DM particles may be produced at 
LHC, and the next-generation direct search experiments will probe a significant portion 



thorough 
If no DM 



of parameter space predicted by various theoretical DM models. In spite o 
DM searches using widely different techniques, the results are negative so far[| 
signature is found in the future experiments, it may suggest that the basic assumption 
that the DM is made of unknown particles is simply wrong. 



x One of the exceptions is the DAMA experiment [3]. However, it is still controversial concerning the 
interpretation of the experimental data. Very recently, two DM-like events were found by the CDMS II 
experiment [I], but more data is clearly necessary to draw definite conclusions. 
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There actually is a DM candidate in the framework of SM, namely, a PBH [5]. In the 
early Universe PBHs can form when the density perturbation becomes large, and it has 
been known that a PBH of mass greater than 10 15 g survives the Hawking evaporation [6J 
and therefore contributes to the DM density [?]. 

In consideration of the entropy of the universe it was pointed out in Ref. [8] that if all 
DM were in the form of 1O 5 M black holes it would contribute a thousand times more 
entropy than the supermassive black holes at galactic centers and hence be a statistically 
favored configuration. Here we consider primordial black holes (PBHs) with masses from 
1O 5 M to 1O~ 8 M and, subject to observational constraints, any of these masses can 
comprise all DM although the entropy argument favors the heaviest 1O 5 M mass. 

There are several ways to realize large density fluctuations leading to PBH formation. 
One is phase transition involving violent processes like bubble collision [9j [TOj [TT] or the 
collapse of string loops (T2j [131 El EE]. As we will see in the next section, however, 
both scenarios have difficulties. Another possibility is the production of PBHs from 
density fluctuations generated during inflation. Since the blue spectrum with a spectral 
index n s > 1 is disfavored by the WMAP data [lj, a single inflation may not be able 
to produce large density fluctuations at small scales unless some dynamics is introduced 
during inflation. On the other hand, the density fluctuations can be easily enhanced at 
small scales in a double inflation model, as first discussed in Ref. [H] in the context of 
the PBH formation. 

In this letter we discuss a double inflation model that consists of a smooth-hybrid 
inflation [T7] and a new inflation [18]. The smooth- hybrid new double inflation was studied 
in Ref. [19] in the context of explaining the large running spectral index suggested by the 
WMAP 1st year data [2D]. In this set-up PBHs with a narrow mass distribution are 
formed as a result of an explosive particle production between the two inflations [21] . We 
will show that the PBH mass can take a wide range of values from 1O _8 M up to 1O 5 M . 
Also, the resultant PBH mass has a correlation with running of spectral index, which 
was roughly estimated in a semi-analytical method in Ref. [22J. Here we numerically 
calculated the correlation, which can be tested by future observations. 

The rest of the letter is organized as follows. In Sec. [2] we briefly review the PBH 
formation and evaporation, and the cosmological constraints on the PBH abundance. In 
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Sec. El we discuss a realistic double inflation model, namely the smooth-hybrid new infla- 
tion model, and show that PBHs with masses satisfying current observational constraints 
are produced. In Sec. 0] we discuss related issues and give conclusions. 

2 PBH formation and observational constraints 

The black hole mass and the formation epoch are related to each other due to the causality. 
In the early Universe, the mass contained in the Hubble horizon sets an upper bound on 
the PBH mass formed at that time. Assuming that the whole mass in the horizon is 
absorbed into one black hole, we obtain 



where Mbh is the black hole mass, M P ~ 2.4 x 10 18 GeV is the reduced Planck mass, 
M Q ~ 2 x 10 33 g is the solar mass, g* counts the light degrees of freedom in thermal equi- 
librium, Pf, Tf and kf are the energy density, the plasma temperature and the comoving 
wavenumber corresponding to the Hubble horizon at the formation, respectively. The 
radiation domination was assumed in the second equality. 

As is well known, Hawking made a striking prediction about the evaporation of black 
holes; any black holes have a temperature inversely proportional to its mass and evaporates 
in a finite time tbh [S], 



Thus the black holes with mass less than 10 15 g must have evaporated by now. PBHs 
which remain as (a part of) DM must therefore be created at a temperature below 10 9 
GeV. In the following we assume that PBHs account for all DM in our Universe. 

The cosmological effects of PBHs have been extensively studied so far. While PBHs 
with masses below 10 15 g are significantly constrained, it is very difficult to detect PBHs 
heavier than 10 15 g because of negligible amount of the Hawking radiation. The MA- 
CHO [23] and EROS [21] collaborations monitored millions of stars in the Magellanic 






(3) 
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Clouds to search for microlensing events caused by MAssive Compact Objects (MACHOs) 
passing near the line of sight. The MACHO collaboration [25] excluded the objects in 
the mass range 0.3M Q to 3OM , and the latest result of the EROS-1 and EROS-2 [26] 
excluded the mass range 0.6 x 1O _7 M < M < 15M , as the bulk component of the 
galactic DM. On the other hand, if we assume that the PBH formation occurs before 
the big bang nucleosynthesis (BBN) epoch, the PBH mass should be lighter than 1O 5 M 
(see Eq. (T5])). Therefore we consider PBHs with masses (i) M BH < 1(T 7 M and (ii) 
3OM < M BH < 1O 5 M . 

Let us comment on other existing constraints. The PBHs with masses heavier than 
43M were claimed to be excluded by the presence of wide binaries [27], but the question 
on the validity of the data used to set the limit was raised by Ref. [28]. Taking account 
of low averaged DM density experienced by the four binaries used in their analysis, the 
strong constraints set by the wide binaries were undermined. Recently, Ricotti, Ostriker 
and Mack investigated the effect of non-evaporating PBHs on the cosmic microwave back- 
ground (CMB) spectrum and anisotropy and found that the PBHs with mass greater than 
~ 0.1 M Q cannot account for the bulk component of DM [2S]- However, the authors made 
assumptions about accretion efficiency in obtaining strong limits on PBH abundances; if 
these assumptions are weakened, all DM could be PBHs for the masses we consider 0. 

The above observational constraints provide us with information on the PBH forma- 
tion. If PBHs are produced at different times, the mass function tends to be broad, 
thereby making it difficult to be consistent with observations. In order to realize the 
PBH mass function with a sharp peak, most of the PBHs should be produced at the 
same time. Thus the production mechanism must involve such a dynamics that only the 
density fluctuation of a certain wavelength rapidly grows. 

What kind of dynamics can create PBHs? First of all, density perturbation must 
become large for PBHs to be formed. There are several ways to realize large density 
fluctuations leading to the PBH formation. One is the phase transition which leads to 
violent processes like bubble collision [91 [101 [11] or the collapse of string loops [HI [131 EH 

2 As pointed out in Ref. [29], the Bondi solution becomes invalid for the PBH of mass Mbh ~ 10 4 M Q . 
In particular, the duty cycle is not well understood because of the complicated feedback effect. If the 
duty cyle is very small, the constraint of |29j can be weakened. 
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[T5] . However, in the case of the bubble collision, the bubble formation rate must be tuned 
to produce the PBH, and the PBH produced from the strong loops tends to have a broad 
mass function. Another possibility is the production of PBHs from density fluctuations 
generated during inflation. In the standard picture of inflation, the inflation driven by 
a slow-rolling scalar field lasts for more than about 60 e-foldings to solve theoretical 
problems of the big bang cosmology. Then no dynamics for producing a sharp peak in 
the density perturbation is expected [30]. However, there is no a priori reason to believe 
that our Universe experienced only one inflationary expansion. Indeed, the cosmological 
gravitino or modulus problem can be relaxed if the energy scale of the last inflation is 
rather low, and it is then quite likely that there was another inflation before the last 
one. If the multiple inflation is a common phenomenon, we expect that explosive particle 
production between the successive inflation periods may produce a sharp peak in the 
density perturbation at the desired scales, which leads to the PBH formation at a later 
time. In the next section, we show that this is actually feasible using a concrete double 
inflation model. 

3 PBHs from preheating 

In this section we provide a double inflation model producing PBHs with a sharp mass 
function as an existing proof. The double inflation model [19j |21j [22] we adopt consists 
of two stages of inflation; the first inflation is realized by smooth hybrid inflation and 
the second one by new inflation. As shown below, the cosmologically relevant density 
fluctuations are generated during smooth hybrid inflation. After the first inflation, the 
inflaton and waterfall fields of the smooth hybrid inflation start to oscillate and decay into 
their quanta via self-coupling and mutual coupling of the two fields. The interesting point 
is that the decays of the scalar fields are largely enhanced through parametric resonance 
and hence the fluctuations of the scalar fields exponentially grow. This process is called 
preheating. During the preheating phase, only the fluctuations at a specific wavenumber 
corresponding to the inflaton mass rapidly grow, and those fluctuations finally turn into 
density fluctuations leading to the production of PBHs with a sharp mass function. The 
role of the second inflation is to stretch the density fluctuations generated during the 
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first inflation and the subsequent preheating phase to cosmologically large scales. In 
the following, we briefly describe the smooth-hybrid new inflation model [for details, see 

Refs.pHIST]]. 

The first inflation is realized by smooth hybrid inflation [31]. In Refs. [HJ [21], the 
smooth hybrid inflation model is built in framework of supergravity and the superpotential 
and Kahler potential are given by 

Wh = S^ + ^^j (m = 2,3,...), (4) 

K H = \s\ 2 + \m\ 2 + \^\\ (5) 

where S is the inflaton superfield, \1/ and ^ are waterfall superfields, fi is the inflation 
scale and M is the cut-off scale which controls the nonrenormalizable term. From the 
above superpotential and Kahler potential together with phase redefinition and the D- 
flat condition, we obtain the scalar potential as 

where a = \/2ReS and ifi = 2Re^f = 2Re^f. Here and in what follows we use the Planck 
unit Mp = 1 and take m = 2 for simplicity. Although the scalar potential fl6]) is derived in 
the framework of supergravity, one may start with ([6]) without assuming supersymmetry. 
The potential (EJ) has a true vacuum at a = and ip = 2i//iM. For a ^ V fiM /2, however, 
the potential for ip has a a-dependent minimum at 

</> ■ ~ (7) 

Note that ip quickly settles down at the minimum during inflation since its mass is larger 
than the Hubble parameter. Then we can integrated out ip and obtain the effective 
potential for a as 

where = v / 2/3 3 ^ 8 (/iM) 1/ ' 4 . If the scalar potential is dominated by the first term, the 
inflaton o slow rolls and therefore inflation occurs. 
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According to the WMAP 5yr data |lj , the curvature perturbation 7Z, the spectral 
index n s and its running dn s /d\iak at the pivot scale fc* = 0.002Mpc _1 are 
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1.031 ±0.055, 
-0.037 ±0.028. 
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From the effective potential (JHJ) we obtain 
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where cr* is the field value of the inflaton when the fluctuation corresponding to the pivot 
scale exits the Hubble horizon. 

The fluctuation corresponding to the pivot scale fc* exits the horizon at t — t* when 
K/a^t*) = H H = fi 2 /\/3 (H H : hubble during the smooth hybrid inflation). Thus the 
scale factor a* = a(t*) is given by 



In a* 



-21n/i - 136. 



(15) 



The e-folding number between the horizon exit of the pivot scale and the end of the 
smooth hybrid inflation is estimated as 



NJar) 



da — 
V 



3a 8 , 



1 



a: 



O* > <T d ) 

(a* < a d ) 



(16) 



where cr e (<C a d ) denotes the field value when the smooth hybrid inflation ends. 

After the smooth hybrid inflation, a and if) oscillate about their minima and decay into 
the a and ip quanta via self-couplings and mutual coupling of the two fields. Since their 
effective masses depend on the field amplitudes and therefore time-dependent, specific 
modes of the a and ip quanta are strongly amplified by parametric resonance. To see this, 
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let us write down the evolution equation for the Fourier modes of fluctuations <7& from 
(ED as 



a'l + 3Ha' k + 



+ + 3m^ t== cos(m a t) 



o k 0, (17) 



a 2 a a ^JlM 

where m a = J8/j?/M and if) is the amplitude of the if) oscillations, (if) ~ \J \iM at 
the beginning of the oscillations.) Neglecting the cosmic expansion, Eq. (jT7J) has a form 
similar to the Mathieu equation which is known to have a exponentially growing solution. 
The detailed numerical simulation showed that the wave number for the fastest growing 
mode is given by 



~ 0.3 m a . (18) 

^osc 

The width of the peak is also determined by the instability band of Eq. (1TT1) . and it is of 
0(O.1)tjv 

The fluctuations amplified by the parametric resonance eventually produce PBHs when 
they reenter the horizon after inflation. The mass of the PBH is approximately given by 
the horizon mass when the fluctuations reenter the horizon. Thus the PBH mass is 
estimated as 

M BH * 1.4xlO 13 M (^^ . (19) 

From Eqs. ffTB"]) and (j!9p the scale factor at the beginning of the oscillation phase is 
estimated as 

lna osc = -114-lnm CT -O.51n(M BH /M ). (20) 
Because the e- folding number N* is equal to lna osc — In a*, we obtain 

JV, = 21 + 0.5 ln(//M) - 0.5 ln(M BH /M ). (21) 

For a fixed black hole mass Mbh, there are two parameters in the model, i.e., \l and M, 
one of which can be removed by using the WMAP normalization ([9]). Therefore observable 
quantities can be expressed in terms of one free parameter, leading to a non-trivial relation 
between n s and dn s /dlnk. In practice, we adopt fiM as the free parameter, and solve 
Eqs. ( IT6l) and (121j) for a* in terms of \iM. Then \i and M are determined with use of 
Eqs. (JI2]) and ([9]) for a fixed \iM. Thus, varying fiM, we obtain sets of model parameters 
which are consistent with the observed curvature perturbations. 



9 



After a and if) decay, the second inflation (= new inflation) starts. As mentioned 
before, the role of the new inflation is to stretch the fluctuations produced during the 
smooth hybrid inflation and subsequent preheating phase to appropriate cosmological 
scales. The effective potential for the new inflation is given by 

Kew = V* (l - °-^) ~ 9 ~V^ + f^ 8 , (22) 

where <f> is the inflaton of the new inflation, v is the scale of the new inflation and g and 
c are constants. The scale factor a/ at the end of the new inflation is estimated as 

where Tr is the reheating temperature after the new inflation. Therefore, the new inflation 
should provide the total e-fold number ~ (In a/ — lna osc ). 

In order to estimate precisely the fluctuations generated in the present model, we 
numerically integrate the evolution equations for the homogeneous modes of a, if), and (f> 
and the Fourier modes of their fluctuations as well as the metric perturbations in the same 
method described in Ref. [21]. In Fig. [T]the power spectrum of the curvature perturbation 
P n is shown for \i = Av = 7.14 x 10" 4 , M = 0.29, c = 0.1 and g = 2 x 10" 5 . It is seen that 
the spectrum has a very sharp peak at k ~ 10 10 Mpc" 1 which corresponds to the PBH mass 
~ 10~ 7 M Q . Since the peak is so sharp, PBHs with very narrow mass range are formed. 
The PBH abundance is determined by the height of the peak which depends on the decay 
rate of the inflaton. We can obtain an appropriate PBH abundance by tuning the decay 
rate [21]. We also predict the spectral index and its running for Mbh = 1O 5 M , 10 _7 M Q 
and 1O _8 M in Fig. |5J Note that the e-folding number of the smooth hybrid inflation is 
smaller than that needed to solve the horizon and flatness problems by a single inflation 
model. Therefore the observable density fluctuations are produced near the end of the 
inflation when the slow roll parameters are larger, leading to a larger value of dn s /d\ak. 

4 Discussion and Conclusions 

In contrast to the conventional WIMP DM model, PBHs have only gravitational interac- 
tions. In order to detect PBHs, we need to carefully look at the effect induced by PBHs 
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Figure 1: Power spectrum of the curvature perturbation. We take /i = 4v = 7.14 x 10~ 4 , 
M = 0.29, c = 0.1 and g = 2 x 10~ 5 . 

such as gravitational lensing, gravity waves, etc. Intermediate mass black holes in the 
range 3OM < Mbh < 1O 5 M can be sought, for example, by higher-longevity microlens- 
ing events [231 121] an d by higher-statistics analysis of wide binaries [271 EH]- in particular 
there appeared recently an interesting idea that if the DM is explained by the PBH of 
mass 10 5 M Q , it may account for the size evolution of the elliptic galaxies by dynamical 
friction [32]. Further observations and theoretical study may reveal the presence of the 
PBH DM. 

What makes the PBH particularly attractive as a DM candidate is that it is naturally 
long-lived due to the gravitationally suppressed evaporation rate. No discrete symmetries 
need to be introduced in an ad hoc manner. Also the PBH DM may be motivated from 
the arguments based on entropy of the Universe [H]. 

In this letter we have argued that the PBH is a natural and unique candidate for 
the DM in the minimal theoretical framework, namely, the SM. Using the smooth-hybrid 
new double inflation model, we have shown that it is possible to produce PBHs of mass 
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Figure 2: n s and dn s /d\nk for PBH mass M BH = 1O 5 M , 1(T 7 M , and 1O- 8 M . 

ranging from 1O" 8 M to 1O 5 M . Importantly, the PBH mass relates the scalar spectral 
index and the running of the spectral index, which can be tested by the Planck satellite. 
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